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TOPICS 

- Summary of paper on spherical harmonics expansion of the cor- 
relation tensor: Rubinstein, Kurien, Cambon; submitted to JoT. 

- Reformulation of the Launder-Reece-Rodi model. 
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SO(3): rotation group in 3D — action on polynomials 



{x, y, z} 

3D: spin 1 




xy,---} = {x 2 + y 2 + z 2 } + {x 2 - y 2 , y 2 - z 2 ,---} 

' " V ' " V ' 

6D id: spin 0 5D: spin 2 

O OOO o 0 2 O 

X y,' • ■} = (x +y +Z ){(x, y, z)} + {(x — xy , x - xz , 

V ' ' V ' ' V 

10D 3D: spin 1 7D: spin 3 



3 = ^2x3 = 3 10 = ±4x5 = 3 + 7 

6 = |3x4=l + 5 15 = j5x6 = l+ 5 + 9 

Single point moments may oversimplify the description of anisotropy. 
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Spherical harmonics basis of 5D space {x 2 — y 2 ,y 2 — z 2 , z 2 — x 2 ■■ ■} 

( x + iy ) 2 = Y 2,2 + iY 2, ~ 2 ( x + iy)z — Y 21 + iY 2, ~ ^ 

3z 2 - (x 2 + y 2 + z 2 ) = y 2 ’ 0 = y 2 


Spherical harmonics basis of 7D space {x 3 — 2xy 2 ,x 3 — 2 xz 2 , • • • }: 

(x + iy) 3 = y 3 ’ 3 + 1 Y 3 - 3 (x + iy) 2 z = Y 3 ’ 2 + iY 3 ~ 2 
( x + iy)(5z 2 — 2(x 2 + y 2 + z 2 )) = y 3,1 + iy 3 ’ -1 
5z 3 - 3 z(x 2 + y 2 + z 2 ) = y 3 ’ 0 = y 3 

SO(3) alone does not provide any basis for these spaces. 
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Correlation tensor Uij(k ) is solenoidal: k l U l .j(k) = k,jU t: j(k) — 0. 
For isotropy, {/y(fe ) = [/(fc)Pj J (fe) where P i; j{k) — 5ij — k~ 2 kikj 
Directional-polarization decomposition (Cambon): 

tfii(fc) = Ufj r (k) + Uff(k) 

where uf- r (k) = ^(U(fe) : P(fe))Pij(fe) and 
[/^(fe) = £^-(fc) - U?j r (k). Or 

= u^ypjjjk) + uff(k) 

tensorially isotropic trace-free solenoidal 
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For axial symmetry (Batchelor, Chandrasekhar) 

Uij(k) = U dir (k)Pij(k) + U pol (k)Sf° l (k) 
where S p ° l (k ) = 

k^cLidj — (a • k)(kidj + kja^ ) + ^(a • fc ) 2 

Spherical harmonics expansions: 

U dir {k ) = ^ A v (k)k~ v Y v (k) 

v> o, even 

y°(fe) = 1, Y 2 (fc) = k 2 — 3 (a • fc) 2 , Y 4 (fc) = 3/c 4 — 30/c 2 (a • fc) 2 + 35(a • k ) 4 

with 

U pol {k) = B u (k)k~ u Z u (k) 

u>2 even 

Z 2 (k) = 7 (a • fc) 2 — k 2 y 4 (fc) = 33(a • fc) 4 — 18/c 2 (a • fc) 2 + k A 


&ij + k 


kjkj 


— ^a 2 k 2 Pij{k) 
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general anisotropy: even spin polarization 


Uff(k) = J2 E A^{k)Y^{k) 

v>0 GVGfl — 

Angular dependence parametrized by amplitudes A u u (k) 

Lj ">r" 

depend only on wavenumber. 

= C'j 3 [Y v ' ,J (k)} 

for rotation invariant (spin 0) matrices of differential operators C\- 
(Arad et al., Zemach, ....) 
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Y^(k) + iY% M (fc) 

= m(m - 1 )(kx + ikyY 


-2 


dkS. 


Y v ik))(Y^ 2 (k)+iY^ 2 (fe)) 


+2/x(fca; + i^)^ 1 (V^O) + iY$ 1 (k)) 

\dk z ) 

+ (k x + ikyY (^2 y "( fe )) ^’ 0 ( fe ) 
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odd spin polarization 

- Rotational strains couple even and odd spins. 

- Similar formulas, but different basis tensors X 2 ’^(/c). 

- Kassinos et al. structure tensor formalism: stropholysis is related 
to spin 3 polarization. The ‘stropholysis spectrum’ 

Qijiik) = e ipq jdS(k ) Uj° l (k)k- 2 k p k e 

defines the projection of U onto its spin 3 component. 
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Reformulation of the LRR model 


Mean flow couplings: 


Uij(k) = - 


U iv (k)^ + U~Jk) dUi 


dx 


p 


dx 


p J 


+ k 


d 

m~ Uij (fc) 


da 


rn 


dk 


n 


dx 


n 


production 


mean— flow distortion 


-|- 2k kikmUpjik) kjkmU v i(k) 


3 


pi 


dU 


m 


dx 


p 


rapid pressure-strain 

Single-point reduction: 

c dJa ' c did 

dk Uij(k) = — dk U iv {k ) — - + / dk 2k~^kik rn U v Ak') 


^ dx 


p 


OK, ‘closed' 


P3 

closure problem 


dx 


P 


Directional-polarization decomposition 


U dir (k) = -U dir (k)(P(k) : S) - U po \k) : S + k m —U dir (k) 

dk n 


dU m 

dx n 



Two models for the correlation tensor consistent with LRR 


Uij(k) = \u{k)P tJ {k) + R^k) 


where, for model I, 


Rij(k) = fe _2 (H(fe) : fcfe)Pj J (fe) 

and for model II, 


k krn^kiHjnj (/c) kjH rn j(k')') 

+ifc- 2 (<5 ii + k- 2 kikj)(H(k) : kk) 


The Reynolds stress deviator is 



dk Rij(k ). 


model I 


Dir-pol 

U dir (k) = U(k ) + k~ 2 V\(k) : kk 
Uf° l (k) = 0 

Stress deviator 

Rij = i 5 J k 2 dk Hij(k) 


Tensor H is basically Kassinos et al. dimensionality. 


LRR and model I 


Mean-flow coupling, directional anisotropy: 


U(k) -)- fc _2 H(fc) : fefe = U(k)k ~ 2 ( S : kk) + fc _4 (H(fe) : fefc)( S : fefe) 


-\-k 


d dUm , , 9 7 _ 2 / 7 \ . , ^dUm 


rrv 


U (fe) — h km— 


dk n 

Evaluating the derivatives: 


k~ z {Y\{k) : fefe) 


<9x 
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£/(fe) + fe _2 H(fc) : fefe = (t/(fe) + fct/'(fc))fe“ 2 S : fefe 

— fe~ 4 (H(fc) : fefe)(S : fefe) + AT 3 (H'(fc) : fefe)(S : fefe) 

+fe“ 2 (H(fe) • S + S • H(fc)) : fefe + fc“ 2 (H(fe) • ft - ft • H(fc)) : fefe 


Spin decomposition: 


U{k) + fc~ 2 H(fc) : kk = (u(k) + kU'{k))k~ 2 S : kk 
spin o spin 2 ' spirit " 

-AT 4 (H(fc) : kk)(S : kk) + AT 3 (H'(fc) : kk)(S : kk) 

S „ ^ S V * * 

spins 2,3,4 spins 2,3,4 

+ /c — 2 (H (/c) • S + S • H(/c)) : fefe + fc“ 2 (H(fc) -fi-Q- H(fc)) : fefe 

' V ' " v ' 

spins 0,2 spin 2 


where 


(H : kk )( S : fcfc) = ||A || 4 + ||A || 2 + ||A || 0 


A || 4 = (H : kk)(S : kk) 


A 

A 


2 = jk 2 {\\ • S + S • H 
0 - ^k 4 ( H : S) 


ffc 2 (H • S + S • H) : kk + ^k A H : S 
§(H : S)l)) : kk 
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Equate terms of equal spin. For spins 0 and 2, 


U(k) = ^(H : S) + : S) 

HO) = (U(k) + kU'(k)) S + f (H(/c) • S + S • HO) - §(H(fc) : S)l) 
+ffc(H / (A:) • S + S • H'O) - §(H'0) : S)l) - (HO) • S2 - S2 • HO)) 

Energy equation: 

k = /°° dA; fc 2 t/ (A:) = /°° dk (j|A; 2 (H : S) + ^A; 3 (H' : S)) = -R : 
Stress deviator equation: 

R = ^kS - ±(R • S + S • R - §(R : S)l) + (R • S2 - Q ■ R) 

This is an LRR model with a special choice of constants. 


However, the spectral equation is not satisfied: spin 4 in the dir 
equation, and spin 2 in the pol equation are both left over. 


Model I cannot satisfy the mean flow coupling equations exactly. 


It is a good approximation provided 

(H : kk )( S : kk ) - ^k 2 (H ■ S + S • H) : kk + ^k 4 H : S « 0 

and 


(U(k) + k ~ 2 ( H : fcfc)) x 

( Pim(k)Pjn(k ) + P in (k)P jm (k ) - Pij(JP)Prnn(X))S 


rrin 


0 


Application of results to assess accuracy/breakdown of LRR model: 
‘in progress’ (at LANL). 


- Results for model II are similar, but spin decompositions are not 
so simple: decomposition of tensor (S : fcfc)H was done by Zemach. 

- ‘Better’ models: including descriptors of higher spins gives a hi- 
erarchy of differential equations. 

- Although stress evolution is determined by spins of all orders, 
stress itself is determined by spin 2 alone. 

- It suggests replacing evolution equations for higher spins by alge- 
braic relations leading to a normal solution : an approximate solu- 
tion of the stress evolution equations in terms of stress-determining 
quantities only. 


